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Abstract 

We consider the problem of evaluation of the weight enumerator of a binary linear code. We show that 
the exact evaluation is hard for polynomial hierarchy. More exactly, if WE is an oracle answering the 
solution of the evaluation problem then p WE — p Ga P p Also we consider the approximative evaluation 
of the weight enumerator. In the case of approximation with additive accuracy 2 an , a is constant the 
problem is hard in the above sense. We also prove that approximate evaluation at a single point e 71 '^ 4 is 
hard for < a < log 2 |1 + w\ » 0.88. 

The weight enumerator of a linear binary code C can be denned as follows 

«>c(«) = 5> |x| - a) 

Here \x\ is a Hamming weight of binary string x S {0, 1}". According to the definition, the weight enumerator 
is a polynomial in q with integer coefficients. The exact calculation of any coefficient of the weight enumerator 
is a #P-problem. Evaluation of the weight enumerator at any positive integer point is a #P-problem too. 
As for evaluation at arbitrary point the approximate setting seems more appropriate. Using approximate 
evaluation with additive accuracy of 2 _poly (") at different points it is possible to compute all (integer) 
coefficients of the weight enumerator. The vector of coefficients is expressed as a linear transform of the 
vector of values in any set of n+ 1 point. A matrix of this transform is inverse of the Vandermonde matrix. 
The accuracy of 2~ poly (") for values of a polynomial degree n is sufficient to get the coefficients of the 
polynomial with 1/2-accuracy. 

To avoid the difficulties that arise when dealing with arbitrary real (and complex) numbers we will 
substitute the problem of the high-accuracy evaluation of the weight enumerator at arbitrary point by the 
problem of computation the coefficients of the weight enumerator. 

Weight enumerator problem. Given a binary linear code C find out the coefficients of the weight 
enumerator wc (<*>). 

We denote an oracle answering the weight enumerator problem by WE. 

Given coefficients of the weight enumerator of a code it is easy to find out the minimal distance of the 
code. The latter problem is NP-complete [8]. Moreover, it is hard to approximate the minimal distance [1]. 
So, the evaluation of the weight enumerator is NP-hard. Here we give more exact characterization for the 
complexity of this problem. We show that an oracle answering a solution for the evaluation of the weight 
enumerator problem can simulate any oracle in the class GapP [2] . The class GapP is closure of #P under 
subtraction. The formal definition of the GapP uses the notion of counting machine. A counting machine is 
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a non-deterministic Turing machine running in polynomial time and finishing at either accepting or rejecting 
halting state. Given an input word x a counting machine produces a gap gM{x). The gap is the difference 
between the number of accepting computation paths and the number of rejecting computation paths. The 
inclusion #P C GapP and Toda's theorem PH C P# p imply that P Ga P p is hard for the polynomial hierarchy 

Theorem 1. P WE = P Ga P p . 

In our reduction we use the techniques of quantum computations. It was proved by Fenner, Green, Homer 
and Pruim [3] that determining acceptance possibility for a quantum computation is hard for the polynomial 
hierarchy. Here we take the similar approach considering the problem of very tight approximation of a matrix 
element (0|?7|0) for a quantum circuit realizing an unitary operator U. 

It will be shown below that quantum circuits can compute the approximation of weight enumerator 
for some codes with an additive accuracy 2 an , where a is constant. So, it is interesting to characterize 
the complexity of approximate evaluation of the weight enumerator with an additive accuracy of 2 an . We 
consider two problems related to this question. 

Estimation at an arbitrary point. The input of this problem is (C, q), where C is an [n, d\ binary 
linear code and q is a complex number e™ v with ip be a rational number. The output w must satisfy the 
condition 

\w-w c (q)\ <2 a ". (2) 
Estimation at the point uj. The input of this problem is (C,q), where C is an [n,d] binary linear 
code and uj — e™/ 4 . The output w must satisfy the condition 

\w-w c (u)\ < 2 an . (3) 

We denote the oracles answering these problems by WE (a) and WE u (a) respectively. 
We will prove that the WE is Turing reducible to both WE(a) and WE u (a). In the latter case we need 
an additional restriction to the constant a: < a < a — log 2 |1 + ui\ w 0.88. 

Corollary 1. p WE ( Q ) = p Ga P p for any < a < 1. 

Corollary 2. P WE -( Q ) = P Ga P p for < a < a . 

These corollaries show that the approximation with exponentially small accuracy of the weight function 
of general linear code by quantum computations is hardly possible. 

For basic facts about quantum computation we refer to the book of Nielsen and Chuang [6]. We will 
need a very little of coding theory. The basic definitions can be found in the classical book of MacWillams 
and Sloane [7] . For details of gap-defined classes see the paper [2] . 

1 Matrix element problem and proof of Theorem 1 

We choose the basis B = {CNOT, H, T} for quantum circuits, where 

CNOT: \a,b)~\a,a®b); H = -±= (| ; T = (J J^. (4) 

It is well-known that this basis provides universal quantum computation. Solovay-Kitaev theorem guarantees 
the fast and effective approximation of any unitary operator by the operators in this basis. 

Matrix element problem. Given a description of a quantum circuit in the basis B and accuracy 
threshold e find out the approximation u of a matrix element (0|t/|0) with accuracy e: \u — (0\U\0)\ < e. 
Here U stands for an operator realizing by the circuit. 

Let's denote an oracle answering the matrix element problem by ME. 

Theorem 2. P ME = P Ga P p . 

To prove Theorems 1 and 2 we construct three Turing reductions: a GapP oracle to the oracle ME 
oracle; then the ME oracle to the WE oracle; then the WE oracle to a GapP oracle. 
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1.1 Reduction an GapP-oracle to the ME oracle 

Let / be a GapP-function. We assume that / is computed by a counting machine in the normal form and 
the machine makes only nondeterministic moves with branching degree 2. Lemma 4.3 from [2] guarantees 
that we do not loss the generality: the gaps computing by machines in normal form and general ones differ 
by factor 2. 

Reduction of /-oracle to the ME oracle reproduce the argument used in [9]. Let q be the number of 
nondeterministic moves along a computation path. The number does not depend on a choice of computation 
path due to the normal form condition. So, we may consider counting machine described above as a deter- 
ministic Turing machine M supplied with an additional input: a guess string r € {0, l} q . Bits of the guess 
string determine the choices for the nondeterministic moves. This machine produces the same gap gKi{x) as 
the initial one. 

Now we transform the machine M to a Boolean circuit Cm in the basis {©, A} computing 1 for pairs 
(input, guess) finishing at the rejecting state and computing for pairs (input, guess) finishing at the accepting 
state. The size of the circuit Cm is quadratic on running time of the M . For fixed input string x the gap 
produced by M is the difference between the number of guess strings u that give the circuit output and 
the number of guess strings u that give the circuit output 1. 

Hereinafter we will simplify notation and will use + and • instead of ® and A respectively 
Let s be the size of the circuit Cm- By z¥\ 1 < k < s, we denote auxiliary variables of the Cm- Wc 
also assume that the circuit output is given by the value of the variable Each assignment in a circuit 
has the form z^ := a * b where * G {+, •} and a, b are either input or auxiliary variables. This assignment 
corresponds the equation Zj: — z k + a * b — 0. Note that the values of input variables and the guess 

string (x, u) determine the values of all auxiliary variables. So, for each x fixed the number of solutions of 

(7) (7) 

the system of equations Z k ' — 0, f < k < s, zf' — Q equals the number of rejecting computation paths of 
the counting machine M for ( = 1 and the number of accepting paths for £ = 0. 
For a polynomial p over the two-element field F2 we denote by Ap the number 



Ap = £(-!)*»>. (5) 



Now let's introduce polynomials over F 2 : 

s-l 

7(3) 

k=l 

(e{o,i},vE{o,iy,we{o,i}. 

Let us check that 



F<(u,z,v) =5>*4 j) +M4 J) +0, F x = (f +w)FM +w(1 + FW), 

(6) 



AF x (u, z, v, w) = AF°(u, z, v) - AF'(u, z, v) = 2 s g m (x). (7) 

The polynomials F^ are linear in variables Vj. So, if values x,u,z does not satisfy a system Z^ = 0, 

f < k < s, z^ = C then all possible assignments to variables v contribute the zero term to the sum (5) 
for the polynomial F!>. If values x, u, z satisfy the system then F^ = for any assignment of v. Thus any 
accepting computation path of M contributes 2 s to AF®(u, z, v) and any rejecting path contributes 2 s to 
AF^(u, z,v). The first equation in (7) follows immediately from the definition (5). 

From the construction described above it is clear that function x F x is computable in polynomial time. 

Now we relate the value AF X to a value (Q\U(F X )\Q) for some unitary operator acting in the space B® N 
of N qubits where N is the total number of variables in the polynomial F x . 

The operator Sj = A J (— I) (controlled phase shift) corresponds to the monomial xj = Ylj eJ Xj of F x . 
Controlled phase shift is defined as 



' A J (-f )|a;i ...x n )= -\x! ...x n ), xj = 1, 

^ A J (— I )\xi . . . x n ) = \xi...x n ), otherwise. 



(8) 
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Let 

n n 

U(F x ) = l[H[j} [] SjY[H\j] (9) 

i=l J€M(F X ) J=l 

where M(-F]c) is the set of the monomials of F x . 
We have 

(0\U(F x )\0) = 

ill J€M(/) Mi,...,ttjv 

= ^F E (-1)^ ( " ) =2- JV AF X . (10) 

Now, we approximate each factor in (9) {/(/) with accuracy S = 2~ N ~ 2 /M by a quantum circuit in the 
basis B. Here M = poly(n) denotes the number of factors in (9). A matrix element (0|f/|0) of a product U 
of these approximations is very close to (0\U (F x )\0) : 

\(0\U-U(F X )\0)\<2- N - 2 . (11) 

So, we may compute the exact value of AF X if an approximate evaluation of (0|C7|0) with accuracy 2~ N ~ 2 
is given. 

To complete the proof of this reduction we note that each factor in (9) acts on 4 bits at most. So, by the 
Solovay - Kitaev theorem [4, 6] the required efficient approximation U can be constructed in poly(n) time. 

1.2 Reduction the ME oracle to the WE oracle 

Consider a quantum circuit in the basis B that realize an operator U. All computational paths of this circuit 
have the same scale factor \j2 N l 2 where N is the total number of Hadamard gates in the circuit. So, to 
compute the value (0|[/|0) we need to sum up phases along all computational paths that start at |0) and 
finish at |0). Note that an Hadamard gate that applies last to the particular qubit doesn't affect the number 
of |0) to |0) paths. Let d be the number of other Hadamard gates. Then there are 2 d summands in the 
expression of (0|£/|0). Each summand is indexed by a binary vector weFj consisted of qubit values just 
after application of the corresponding Hadamard gate. 

A phase shift along a computational path consists of two factors. The first describes the action of T 
operators. It has a form 

t 

[JV fc(M) , where oj = e i7r / 4 , 
fe=i 

and £k{u) is a F 2 -linear form. The second factor describes phase shifts induced by Hadamard gates. It has 
a form 

d 

(-1) B( * U \ where B(u) =Y J b j {u)u j . 

3=1 

F2-linear form bj in this expression is the value of a qubit just before application of jth Hadamard gate. 
Using equations i a(Bb = i a i b (—l) ab and i = uo 2 the second factor can be rewritten in the form similar to the 
first factor. 

Finally, we get the following expression of the matrix element (0|i7|0): 

(OTO^^Ti E f[^ k{u \ (12) 

u( Zfd k=l 
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where (3k are F2-lincar forms. 

It is clear that up to the factor 1/2 w / 2 this expression equals wc(w) — the value of weight enumerator at 
the point u for some binary linear code C. Coefficients of linear forms (3k{u) are columns of the generator 
matrix for this code C. 

So, we reduce the matrix element problem to the weight enumerator evaluation problem. 

1.3 Reduction the WE oracle to a GapP oracle 

Let C be a linear binary code. It is a ci-dimensional subspace of an n-dimensional space F 2 . By usual 
properties of GapP functions [2] the function E(D(C),q) = wc(q) € GapP, where D(C) is a description of 
the code C and q is integer. The required reduction follows from the fact that all coefficients of wc(q) can 
be restored efficiently from the value E(D(C),2 n ). 

2 An approximate evaluation 

Theorem 1 establishes the hardness of high-accuracy approximation of the weight enumerator. 

Now we will prove two results concerning the approximation with moderate accuracy. Our choice of accu- 
racy bounds is motivated by observation of Knill and Laflammc [5] that the power of quantum computation 
is just the ability to produce a Bernoulli sequences with a parameter p given by a real or imaginary part of 
a matrix element (0|f/|0) of an operator realized by a quantum circuit. So, value of the matrix element can 
be estimated with an accuracy n~ 0( and exponentially small error probability. 

The formula (12) shows that the value of weight enumerator t«c(w) can be evaluated in some cases with 
rather high accuracy n~°^2 N ^ 2 . Interesting quantum computation that does not reduced directly to the 
probabilistic one corresponds the case N < 2d/\ogd. 

Corollaries 1 and 2 show that the general case differs from the case of these special codes. In this section 
we give the proofs of these corollaries. 

2.1 Combining codes: direct sums and wreath sums 

In this section we introduce two operations with linear codes and describe corresponding transformations of 
weight enumerators. 

Consider an [ni,c£i] linear code A and an linear code B. Let G A and Gb are the generator 

matrices of codes A and B respectively. 

A direct sum A © B is a [m + n 2 , d\ + d 2 ] code defined by generator matrix as follows: 

(t I) ■ 

The weight enumerator of the direct sum is just the product of weight enumerators of summands: 

WA®B(q) = w A (q)w B (q). (14) 

A wreath sum A\B is a [n\n 2 , d\ + d 2 ] code. To define elements of the generator matrix Gaw for AlB we 
assume that rows of Gaw are indexed by pairs (fci,0), 1 < fci < d\, and (0, fc 2 ), 1 < fc 2 < e? 2 , and columns 
are indexed by pairs (Zi, I2), 1 < h < n-i, 1 < h < n 2 - In this setting the elements of Gaib are given by 

{GAlB)(kfl){hM) = { G A)k 1 l 1 (15) 

{G A iB)(o,k)(h,i 2 ) = (Gs)fc 2 i 2 . (16) 

Each pair of words (a, 6), a e A, b e B determines in natural way a word c a b € A I B. Each bit a in the 
codeword a is extended to a string a' 2 . This string is added by the codeword b to produce the Z 2 -block of 
the codeword c a ^ h . So, for weights of these codewords the following relation holds 

|c fe| = \a\(n 2 - \b\) + \b\(m - \a\). 
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From this equation we get the following relation between weight enumerators: 



WAw(q) = £ <f l2|a W(l + q ni - 2H ). (17) 

aeA 

2.2 Proof of Corollary 1 

We will construct a Turing reduction of the weight enumerator evaluation problem to the WE(a) problem. 

Namely, we will evaluate weight enumerators of a code C\ — C Ik at several point qj . Here Ik means 
the trivial code given by identity generator matrix. The parameter k will be chosen later. 

Choose r such that j3 — 2 Q /(2 — r) < 1. Note that (2"( d+fc ) Revaluation of wc k implies much stronger 
estimation of wc- if \q — 1| < r and \wc k (q) — w(q)\ < 2 Q ( d+fc ) then 



wc(q)(l + q) k -w 



w 

w c (q) - 



(! + «)» 



< 2 ad (3 k = 2-°< fe ). (18) 



(In asymptotic notations we assume that k — > oo while other parameters are fixed.) 

To compute the coefficients of weight enumerator wc (q) we use this tight bound for n + 1 points inside 
the circle \z — 1| < r. The distances between these points can be chosen greater than r/(2(n + 1)). So, the 
inverse of Vandcrmonde determinant for these points is 2°(" lo s™). 

Thus, for some k = poly(n, d) we can compute the coefficients of the weight enumerator with accuracy 
1/2. 

2.3 Proof of Corollary 2 

Note that reduction in the subsection 1.2 uses the evaluation of the weight enumerator at the point lu only. 

To construct a Turing reduction of the weight enumerator evaluation problem to WE W (a) problem we 
will obtain an evaluation of the weight enumerator of a code C at point u using the (2 Q (" +a+fe ))-cstimations 
of weight enumerators for codes C(a,k) — C a I Ik at the same point. Here Ik stands for a trivial code of 
dimension k and code C a is produced from the code C by padding a zeroes to each codeword. From (17) we 
get for these codes: 

^c(a,k) - £ u; fe H(l + w »+- 2 H)fc. (19) 
toec 

We assume that fc = 1 (mod 8). 

Divide the weight enumerator sum for code C in four parts: 

w c (uj) = M + Mi + M 2 + M 3 ; M,- = £ u;- 5 '. (20) 

w£C, \w\=j (mod 4) 

It is clear that \Mj\ < 2". 
Now rewrite (19) as 



wc(a.k) = M (l + w" +a ) ft + Mi(l - i^" +a ) fc + M 2 (l - uj n+a ) k + M 3 (l + ^" +a ) fc (21) 

Take a such that n + a = 1 (mod 8). Let w\ be an (2 a (" +a+fc ) Revaluation of the Wc( a ,k) an( i ^ = 
From (21) we have 

\M + M^- 1 + M 2 il^ + M 3 (1 ~^ fc " 1)fc - < 2«(«+ a + fc )/iV fc = 2"<"+"> ^) " . (22) 
The last factor is 2-°<' i; ). Note that 

ATfc and N k ~ z 
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So, for some k = poly(n, — loge) we obtain an e-evaluation fio of Mq + M\u)~ . 

Incrementing a by 2 in the above argument we get e-evaluations fij of Mj + Mj+iu; -1 . Solving the linear 
system 

x 1+ x 2 ^=^ 

X2+X3LO =yU 2 , 
X 3 + XqUJ- 1 = /i 3 , 

we get the 0(e)-approximations for Mj. It is enough to solve the weight enumerator evaluation problem for 
code C. 
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